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Polarization dependence of emission spectra of multiexcitons in self-assembled 

quantum dots 



N.Y. Hwang and S.-R. Eric Yan£| 
Physics Department, Korea University, Seoul Korea 

We have investigated the polarization dependence of the emission spectra of p-shell multiexcitons 
of a quantum dot when the single particle level spacing is larger than the characteristic energy of the 
Coulomb interactions. We find that there are many degenerate multiexciton states. The emission 
intensities depend on the number of degenerate initial and final states of the optical transitions. 
However, unlike the transition energies, they are essentially independent of the strength of the 
Coulomb interactions. In the presence of electron-hole symmetry the independence is exact. 
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I. INTRODUCTION 

During recent years great progress has been made 
in the development of experimental tools for study 
of optical properties of single dot structures Con- 
focal microscopy and near field scanning microscopy 
have been particularly useful. Self-assembled quan- 
tum dots(SAQDs) are promising as single photon 
generators 0, Q because of the high optical quality of 
the resulting dot structures. Photons may be generated 
from states formed by a number of electrons and holes, 
i.e., multiexcitons |, E 1 0, & @, El El El El EI 
El El El El El- Fine structures of excitons, biexci- 
tons, and triexcitons have been studied experimentally 
0, 0, H 0, B El EH • Exact numerical diagonalization 
and quantum Monte Carlo studies of emission spectra of 
multiexciton states have been performed for several dif- 
ferent shapes of SAQDs[H El El, El, El El El, El- 
The emission spectra of Bose-Einstein condensed magne- 
toexcitons in the strong magnetic field limit have been 
also investigated theoretically in SAQDs[20(. 

However, the polarization dependence of the emission 
spectra of multiexcitons in quantum dots is not well un- 
derstood. In one-, two-, and three-dimensional electron- 
hole plasmas with parabolic bands the ratio between 
the emission intensities for linear and circular polariza- 
tions is constant since their optical strengths are M 2 and 
2M 2 , see Table Q] (The constant M is the dipole ma- 
trix clement — e(x|x|s). where |.r) and \s) are the Bloch 
wave functions). This is because the independent quasi- 
particle picture holds in these systems and each single 
particle state is always doubly spin degenerate. How- 
ever, in zero-dimensional quantum dots this may not be 
the case since single particle levels are discrete and many 
body effects are strong, see Figs. [T]and[2] It is unclear to 
what extent the independent quasi-particle picture holds. 
A simple example demonstrates that there may be some 
polarization dependence. The emission from the s-shells 
is simple to analyse and shows that the emission spectra 



for linear and circular polarizations are not proportional, 
see Fig. [5] This is because for N = 2 — > 1 and linear 
polarization the initial state of the optical transitions is 
non-degenerate while the final states are doubly degen- 
erate. This means that the emission strength doubles 
from M 2 to 2M 2 . We investigate polarization depen- 
dence of the emission spectra of a quantum dot when a 
photon is generated by the recombination of an electron- 
hole pair in the p-shells, see Fig. [TJ In this case the 
polarization dependence of the emission spectra is ex- 
pected to be more complicated since many body states 
can be strongly correlated. Our results demonstrate that 
the optical strength can depend sensitively on the polar- 
ization state of the emitted photon. The optical strength 
can be enhanced significantly if there are several degen- 
erate multiexciton final states. In addition the optical 
strength can take different values depending on which de- 
generate initial state of the optical transition is occupied. 
The emission spectra are calculated using these optical 
strengths, and are given in Figs. [3] and |U The emission 
intensities are almost independent of the strength of the 
Coulomb interactions in our model, although the transi- 
tion energies do depend on the Coulomb interactions. In 
the presence of electron-hole symmetry the independence 
is exact. Our results may be useful since they could pro- 
vide an important qualitative result that must be valid 
in a more general approximation. 
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TABLE I: Transition matrix elements between a heavy hole 
and an electron for different polarizations. Note that a hole 
is the absence of a valence electron, and its spin takes the 
opposite value of the valence electron. 



II. MODEL HAMILTONIAN 

For our purpose we may take a simple model, which is 
illustrated in Fig. [T] In lens-shaped self-assembled quan- 
tum dots the conduction and valence band electrons can 
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FIG. 1: Emission from the p-shells when N pairs of electrons 
and heavy holes occupy s- and p- orbitals of two-dimensional 
harmonic oscillators. Note that both electron and hole en- 
ergies are measured positive. The bandgap Eg defined as 
the energy difference between the bottom of the electron har- 
monic potential and the top of the hole harmonic potential. 
The numbers 1 and —1 stand for the angular momentum com- 
ponents along the axis perpendicular to the two-dimensional 
layers. 
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e^. Electron wavefunctions are u^rjx^^) f fX^) and 
heavy hole wavefunctions are u v ,±i(r)x(±i)g£(f). Here 
the conduction and heavy hole valence band Bloch wave- 
functions are denoted by u c (r) and u v ^±\(r). In the fol- 
lowing we will represent hole states u v ^i{r)x{\)gi{'r) and 
M Vj _i(r)x(— \)ge.{r) as having spin components 1/2 and 
— 1/2. The two-dimensional electron and hole effective 
mass wavefunctions are fe{fi) and gi(?2)- The quantities 
and at = „ » n arc the characteristic length 

n 2m: ii. o 



2m* Sl c "h 2mjf2 h 

scales for electrons and holes. The parameters m* and 
m£ are the effective masses of electrons and holes, and 
f) e and fih are the strengths of the harmonic potentials. 
Two-dimensional polar coordinates are n, T2 and 8. We 
assume that the wavefunctions have a width d along the 
axis perpendicular to the two-dimensional plane. This ef- 
fect is included in the shape of the Coulomb interactions. 
The quasi-two-dimensional Coulomb interaction has the 
form V{r\,r2) = — -. e , where d is the width of 

e V \r 1 -r 2 \ 2 +d 2 

the single particle wavefunctions along the z-axis (e is 
the background dielectric constant). Single particle en- 
ergies of electrons and holes are e\ — M7 C (|^| + 1) and 



1). Note the hole energy is defined to be 



positive. The bandgap Eq, defined as the energy differ- 
ence between the bottom of the electron harmonic po- 
tential and the top of the hole harmonic potential, is set 
to zero. In our model we assume that the single particle 
level spacing is larger than or comparable to the char- 
acteristic Coulomb energy, which is often valid in lens- 
shaped SAQDs[21]. This allows us to ignore excitations 
from s to p levels and from p t o d levels/continuum states 
of two-dimensional layers [2l|, [22[ • We can thus assume 
that the s-shell is completely filled and inert. 
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FIG. 2: Optical transitions between the s-shells. Filled (open) 
triangles represent electron (hole) spins, (a) Possible transi- 
tions for circular (P + ) and linear (P x ) polarizations when 
N = 1 — * 0. (b) Possible transitions when N = 2 — * 1. (c) 
Schematic display of the emission spectra for right circular 
polarization, (d) Emission spectra for linear polarization. 



be approximated by those of harmonic oscillator wave- 
functions. Since in quantum dots only states near k = 
are relevant we will neglect the heavy and light hole mix- 
ing and concentrate only on the physics of the heavy 
hole. The typical energy spacing between single particle 
energies is 20 — 40meV, which is larger than the strength 
of the Coulomb interaction of order lOmeV. The ener- 
gies of s-level and p-level orbitals are denoted by e^' h and 



The many body Hamiltonian H for the electron- 
hole system can be written as a sum of single par- 
ticle terms Hq and interaction terms VS n t, i.e., H = 
H a + V int with H a = H c + H h and V int = H cc + 
i?hh + -ffch, where the single particle Hamiltonians 
H c = J2e,a e e a lcr a ^ and H h = J2e,a e ^L^- Tlic 
operator a\ and b\ create an electron and a heavy 
hole (o =t or J,). Electron-electron, hole-hole, and 
electron- hole interaction terms are H cc , ii/hh, and -ff e h, 
respectively. The electron Coulomb matrix elements 

„ r „ TTCC 



5ai ,o-4 3<r 2 ,<r 3 titi+ti M+U (ftife 2 \y\ fii fe 3 ) 



and hole Coulomb matrix elements are f^gJ/3 2J 3 4/ 3 3 



0~(7\ ,CJ4 3(72 ,^3 



s e 1 +£ 2 .i a +e i (ge 1 g£ 2 \V\gi i gi 3 } (a and f3 stand 
for (I. a)). The electron-hole matrix elements are 
U fh<x 2 f3 4 a 3 = ^ 6 TwrJcr 2 ,<T 3 Se 1 +e 2 4 3 +e i (gi 1 fe 2 \V\ge i f e;i ). 
Since the single particle energies are independent of spin 
we can write e\ a = e\. When S ai _ (Ti = S a2iCF;s = 1 we may 
write Ut?< ri J a n,u< Ul t a e a = ^AA' Similarly, simpler 
notations may be introduced for U hh and U eh . 
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III. MULTIEXCITON HAMILTONIAN MATRIX 

The total z-component of angular momentum L z = 
L c .z + -^h.z is conserved. The electron and hole angular 
momenta, L oz and Lh,z, need not be conserved sepa- 
rately since H h can change them. On the other hand 
the total z-component of spin for electrons (holes), S e ^ z 
(<Sh,z), is conserved since i? h cannot change S c , z or Sh, z - 
If S c ^ z is known then the electron spin quantum number 
S e can be deduced: When S c . z = —1,0, or 1 then S e = l, 
and when S c . z = -| or | then S c = \. In our model 
S'o.z cannot take other values than the ones just listed, 
see Fig. [5] Similar results hold for hole spins. We label 
the Hubert subspaces by \L Z , S e ,zi Sh, z ). The resulting 
Hilbert subspaces are listed in the Fig. [5] (see Appendix 
B). 

The basis vectors in each Hilbert subspace may be 
chosen as the single Slater determinant states \<fii) — 

a l N ■ ■ - a k b l N ' ■ ' 6 kl )' with i = ("i'A) where a, = 
(ai, 0C2, aisr) and — (/3i,/?2, 0n)- The kth eigen- 
state of N electron- hole system is written as |$fe) = 
J2i Ai\<f>i). The eigenstate |$fe) and eigenvalue Ef. satisfy 
the following matrix equation 'Ylij{4>i\H\4>j)A^ = E^A\. 
The diagonal elements of the Hamiltonian matrix are 

(<l> i \H\ ( l> i )=E e + E h +E eh . (1) 

The total energy of the electrons/holes can be determined 
from 



(2) 



a/3 



with L7*j = Vj* af) - U^ ppa , where p = e,h. The elec- 
tron and hole occupation numbers in |0j) are denoted by 
n e a and n\. The total electron-hole interaction energy is 
Ed, = E„(3 n a n ^5a^ Tne off-diagonal elements are 



M#l0 3 ) 



\Hee\4>3 



|#hh|0j 



i\H e h\(f)j) 



(3) 



IV. EIGENSTATES 



In our model it is possible to find analytically all the 
eigenstates in each Hilbert subspace. This is because 
the maximum number of basis states of these Hilbert 
subspaces is at the most six and the Hamiltonian matrices 
have simple forms. 



A. Uncorrelated single Slater determinant states 
and correlated states with two basis vectors 



In our model many Hilbert subspaces are one- 
dimensional. They are displayed in Fig. [5] and their 



energies are given by Eq. ([T]). There are many two- 
dimensional Hilbert subspaces. Their Hamiltonian ma- 
trices are symmetric and the diagonal elements are the 
same: 



(4) 



(5) 



(falHifa) (falHlfc) \_f-y s 
{h\H\4>x) (h\H\h) ) \5 1 

The eigenvalues of this matrix are 

Ei = 7 + 5, E 2 = 7 - 5, 
and the corresponding eigenvectors are 

|*i) = -±=(|0i) + |0 2 )), |$ 2 ) = - |0 2 )). (6) 

Note that these eigenvectors are independent of the 
Coulomb interactions. The off-diagonal element 5 — 
(0i \H \cf>2) is non-zero for N = 3,4,5. 



B. Correlated states with six basis vectors 

In our model there is only one Hilbert subspace with 
six basis states: the Hilbert subspace \L Z , S C:Z , Sh, z ) — 
1 0,0,0) for N — 4. The basis vectors are de- 
fined as follows: |0i) = a\ ^a_ x ^b\ ^b'^ AS), \<j) 2 ) — 

t „t ' J i,t ' 



4^ ^ ^ 

4^ a -i,ri,i b -i,i\ s )< I&) = a -i,i a ~i^ b l,i b l,-\\ s )' I0e) 

a\ <4 t b -i L-i l\S)> where \S) = 4 ^4 ^ 6q ^b ||0) rep- 
resents the filled s-shells. Note that | 0b) and |0e) con- 
tain doubly occupied p-orbitals. The diagonal matrix el- 
ements are A = (4>k\H\4>k) and are given by Eq. |T]). 
If we define B = -V^ hl _ Xi C = -V^,i,i-i and 
D = V~[ 1 1 _i, the Hamiltonian matrix is given by 
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with the eigenvalues Ei = A + B + C,E 2 = \{2A - 
B~C- y/B 2 +2BC + C 2 + 32D 2 ), E 3 = A + B - C, 
E 4 = A, E 5 = A - B + C and E 6 = \(2A - B -C + 



Vb 2 



$2) 
$3) 

$4) 
$5) 

*6> 



2BC + C 2 + 32D 2 ). The eigenvectors are 

= 5(|0l) + 102) + |03) + |04)), 

= ai|0i) + a 2 \4> 



5H0O 
i(-|0x) 

&l|0l) + ' 



h) + a 3 |03) + a4|' 

h) - |0 3 ) + |0 4 )), 

106)), 

h) + |0 3 ) + 104)), 
>2> +&3|0 3 ) +b 4 \(f. 



4) + a5|0s) + ae|06), 



(8) 



Note again that |$i), l^a), |3>4), and |$s) are independent 
of the Coulomb interactions and model parameters. On 
the other hand , the expansion coefficients of ($2) and 
I $6) do depend on them. The expansion coefficients can 
be obtained analytically but their expressions are rather 
complicated and lengthy. The lowest energy state is 
since B < 0, C < 0. In this state both electrons and 
holes are in a spin triplet state S e = 1 and Sh = 1 with 

{S e ,M,S h ^) = (0,0). 



C. Number of groundstate degeneracy 

Our investigation shows groundstates have L z = 0, see 
Fig. O When TV = 6, 5, 4, 3, 2 we see from Fig. [S]that the 
number of degenerate states with L z — are 1, 4, 9, 4, 1. 
Let us examine why the degeneracy is 9 for N — 4. From 
symmetry considerations one can predict the number of 
groundstate degeneracy. Confirming this number of de- 
generacy provides a powerful check of the correctness of 
our analytic calculation of eigenenergies. We find that 
the lowest energy states are the ones with electrons and 
holes in triplet states. The electron and hole spin num- 
bers are thus 5*0 = 1 and Sh = 1. This implies that the 
possible values of (S c , z , Sh,*) are 3x3 = 9. Since the in- 
teractions are spin invariant these states are degenerate. 
Our calculation shows that the groundstates of the fol- 
lowing nine Hilbert subspaces |0, 1,1), |0, 0,1), |0, 0, — 1), 
|0,-1,-1), |0,-1,1), |0,1,-1), |0,-1,0), |0,1,0), and 
1 0, 0, 0) are degenerate (see Fig. [5]). Other degeneracies 
may be shown in a similar manner. 



OPTICAL STRENGTH AND EMISSION 
SPECTRA 
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TABLE II: In our model most optical dipole strengths can 
take only four different values, ±M 2 , M 2 , 2M 2 and 4M 2 . The 
oscillator strengths of circular polarization are twice those 
of linear polarization, just like in the case of single particle 
transitions. 

In this paper we consider linearly, right circularly, and 
left circularly polarized photons. We assume fast energy 
relaxation so that only emission from ground states of N 
exciton pairs are needed. The emission spectrum is given 
by 



/(«) = £ K*/|i|*G>l 2 s(£i -E f - u), 
f 



where the luminescence operator is 



L=^2(-t,-a\i-0[£,<r)b. 



(9) 
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FIG. 3: Emission spectra for linear polarization P x . Closely 
spaced vertical lines under the curly bracket have all the same 
transition energy. The numbers above the arrows denote the 
number of different initial states giving rise to the same tran- 
sition energy. The change in the number of electron-hole pairs 
TV — > N — 1 is indicated in the upper part of each panel. 
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FIG. 4: Emission spectra for circular polarization P + . 



An electron and a hole can only recombine when they 
have opposite values of the z-component of angular mo- 
menta and spins. This leads to the selection rules 
AL Z — and AS Z — 0. The optical dipole strength 
from a ground state of N-exciton system, $>q(N), to a 
ground/excited state of (N-l)-exciton system, $/(TV— 1), 
is |($/|L|<I>g)| 2 . The possible values of the optical dipole 
strength are displayed in Table HIl In our model there 
are only three types of dipole transitions: (a) transitions 
between single and two-basis states, (b) transitions be- 
tween two-basis states, and (c) transitions between two- 
and six-basis states. We have computed emission spec- 
tra for polarization along the x-axis (P x ), left circular 
polarization (P~), and right circular polarization (P + ). 
They are shown in Figs. Oand[3J The following param- 
eters are used: ah = 0.46a G , d — 0.8a c , hfl c = 275c, and 
Oh = 0.7f2 G . In our calculation we measure energy with 
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respect to the Coulomb energy Ec = Note that the 
transition energy is actually Ei — Ef + Eg, but we will 
drop the bandgap energy Eq for convenience hereafter. 
The emission spectra for right and left circular polariza- 
tions are identical. Note that only one transition energy 
is possible when N = 6 — > 5, 4 — > 3, and 3^2. This 
is because transitions from the ground state of N = 6 to 
the excited states of N = 5, and from the ground state of 
TV = 4 to the excited states of N — 3 are optically forbid- 
den or vanish. On the other hand, for N = 5 — ► 4 three 
transition energies are possible. One can show that the 
excited states |$3), |<&4) and |$s) of the Hilbert subspace 
|0,0,0) for N = 4 do not contribute to optical strengths 
(see Appendix A). Only |$i), |<&2) and |<&e) contribute. 

Let us discuss the emission spectra of 5 — > 4 for 
P 1 . The oscillator strength of the lowest energy peak 
6.39i?c is 0.0088M 2 . It originates from the transitions 
|0,i,-i)i or |0,-i,i)i -> |$ 6 ), see Fig. [5] (These 
are defined as type A transitions). The subscript 1 in 
|0,|,—5)i indicates that it is the lowest energy eigen- 
state. This oscillator strength is negligible compared to 
those of the others. The oscillator strength of the next 
lowest transition energy 6.58Ec is 1.4912M 2 . It orig- 
inates from the transitions |0, 5,— 5)1 or |0, — -g, 5)1 — > 
I $2) (These are defined as type B transitions). The 
dependence of the strength of this peak on the ratio 
fih/^c is negligible: as it is reduced from 0.7 to 0.1 
the oscillator strength changes smoothly from 1.4912M 
to 1.5000A/ 2 (the change is only about 1%). The de- 
pendence on the Coulomb interactions, i.e., on the ra- 
tio d/a e is also negligible: as it is reduced from 0.8 
to 0.2 the oscillator strength reduces continuously from 
1.4912M 2 to 1.4780M 2 . As these parameters change 
the sum of the two oscillator strengths of transitions 
A and B remains a constant equal to 3/2M 2 . This 
is evidence that our calculations are correct. The 
highest peak with the energy Q.&IEq has the opti- 
cal strength 5/2M 2 and originates from the transitions 
|0,|,-|)i -> |0, l,-l)i and |$i), or |0,-|,|)i 
jO, — 1, l)i and |$i), see Fig. [SJ In these transitions there 
are two possible final states with the same energy, and 
the oscillator strength is the sum of the two contribu- 
tions. The other peak at the same transition energy 
with the strength 2M 2 originates from the transitions 
|0,|,|)i -► |0,l,0)i and |0, 0,l)i, or |0,-i,-|)i -> 
|0, — l,0)i arid 1 0, 0, — 1) 1 • Note that some transition en- 
ergies can have more than one optical strengths. For 
example, for N — 5 — > 4 at the transition energy Q.QIEq 
two values of optical strengths 2M 2 and 5/2M 2 are pos- 
sible. In this case there are two different initial states 
giving rise to the same value of the optical strength 2M 2 . 
The same is also true for 5/2M 2 , see Fig. [31 

Now we discuss the emission spectra for P + for the 
case N = 5 — > 4. The oscillator strength of the lowest 
energy peak 6.39£ c is 0.0176A/ 2 . It originates from the 
transition |0, — A, |)i — > \Q>$), see Fig. [5] (We define it as 
type C transition). This oscillator strength is negligible 
compared to those of the others. The oscillator strength 



of the next lowest energy peak 6.58E C is 2.9824M 2 . It 
originates from the transitions |0, — 5,5)1 — > ($2) (We 
define it as type D transition) . Again as the physical pa- 
rameters change the sum of the two oscillator strengths 
of transitions C and D remains a constant equal to 3Af 2 . 
Note that some transition energies can have more than 
one optical strengths. For example, for iV = 5 — ► 4 
at the transition energy 6.61i?c three values of optical 
strengths M 2 , 2 AT 2 and 4M 2 are possible. The highest 
peak 6.6lEc with the strength 4M 2 originates from the 
transitions |0, 5,— |)i — > 1 0, 1 , — 1) 1 - The peak 6.61E C 
with the strength 2M 2 originates from the transitions 
|0,i,i)i -> |0, l,0)i, or |0,-5,-±)i -> |0, 0,-l)i. The 
peak 6.61Ec with the strength M 2 originates from the 
transitions |0, — 5, 5)1 - * |$i)- There are thus two dif- 
ferent initial states giving rise to the same value of the 
optical strength 2A/ 2 , and there are only one initial states 
giving rise to the optical strengths M 2 and 4M 2 . 

We have verified numerically that when electrons 
and holes have exactly the same properties, i.e., when 
electron-hole symmetry is present, the dependence of the 
oscillator strengths for the transitions A, B, C and D 
on the Coulomb interactions disappears completely. The 
oscillator strengths of these transitions become, respec- 
tively, 0, |, 0, and 3. 

VI. DISCUSSIONS 

Our main results may be attributed to the presence of 
degenerate states and the strongly correlated nature of 
multiexciton states. In the non-interacting model, spin 
properties of electronic states do not change as a func- 
tion of energy. In the interacting case, spin properties 
of electronic states do change as a function of energy. 
At different transition energies spin properties of opti- 
cally active initial and final states are different. Which 
initial and final states are connected by the optical selec- 
tion rules will depend on the polarization of the emitted 
photon. Several factors contribute to the quantitative 
aspects of the polarization dependence of the emission 
spectra. As illustrated in Fig. [2] the presence of degener- 
ate final states can enhance the oscillator strengths. The 
presence of degenerate initial states can lead to the pres- 
ence of degenerate transitions energies. The correlated 
nature of some eigenstates influence the value of their 
oscillator strengths and can lead to a cancellation of the 
optical strength. 

Our approximate calculation shows that there is a 
significant polarization dependence. For more realistic 
SAQDs our results may be used as a starting point of 
a perturbative treatment of the electron-hole exchange 
interaction, the coupling between s and p and between 
p and d levels, coupling between p-lcvcls and the con- 
tinuum states of the two-dimensional layers, heavy and 
light hole coupling, and elliptic dot confinement potential 
|23j |. We expect in these cases that the emission intensi- 
ties will depend on the Coulomb interactions. However, 
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as long as the single particle level spacing is larger than 
the characteristic Coulomb energy the effect would be 
minor. 
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